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Research  under  this  grant  focussed  on  several  aspects  of 
gravity  wave  excitation,  propagation,  and  dissipation  that  are 
expected  to  be  important  in  the  atmosphere.  Initial  studies 
addressed  the  excitation  of  propagating  waves  by  unstable 
shear  layers  and  found  that  the  nonlinear  interaction  of 
evanescent  unstable  inodes  is  an  efficient  source  of  such 
motions.  Other  numerical  studies  examined  the  consequences  of 
gravity  wave  propagation  and  saturation  in  the  middle 
atmosphere.  Important  findings  include  an  amplitude  limit 
imposed  by  wave  field  instabilities,  the  self -acceleration  of 
large -amplitude  motions  which  may  greatly  expand  the  phase 
speed  distribution  of  mesospheric  wave  motions,  and  the 
relative  insignificance  of  nonlinearity  in  limiting  wave 
amplitudes  and  preventing  wave  field  instability. 
Observational  studi*?:  revealed  wave  field  dynamics  to  be 
largely  consistent  with  linear  instability  theory,  with 
turbulence  produced  at  that  s-ite  in  the  wave  field  where  the 
motion  is  most  unstable.  Wave  amplitudes  were  seen  to  be  near 
saturation  values  and  easily  described  by  a  simple  saturation 
model  of  the  evolving  gravity  wave  spectrum  throughout  the 
atmosphere.  Additional  theoretical  studies  addressed  the 
turbulent  transport  of  heat  and  constituents  and  the  induced 
mean  vertical  motions  due  to  vertically  propagating  gravity 
waves,  contributing  to  our  understanding  of  apparent 
differences  between  observations  and  modeling  results. 
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1 .  INTRODUCTION 

Research  supported  under  Grant  AFOSR  82-0125  has 
contributed  substantially  in  an  area  of  emerging  importance  in 
atmospheric  dynamics  over  the  past  four  years.  This  area  is 
the  study  of  atmospheric  gravity  waves,  including  their 
excitation,  propagation,  dissipation,  transports  of  momentum 
and  energy,  generation  of  turbulence,  and  diffusion  of  heat 
and  constituents.  Also  of  importance  are  the  very  significant 
effects  of  such  motions  on  the  large-scale  circulation  and  the 
thermal  and  constituent  structures  of  the  lower  and  middle 
atmosphere.  While  this  field  has  benefitted  from  contributions 
by  many  authors  dating  back  to  the  pioneering  works  by  Hines 
(1960,  1963,  1972),  Hodges  (1967,  1969),  and  Bretherton  (1966, 
1969a,  b),  the  purpose  of  this  report  is  to  document  the 
contributions  of  this  research  effort.  We  will,  therefore, 
only  refer  to  other  contributions  where  they  indicate  the 
relevance  of  the  present  results. 

The  evolving  objectives  of  this  research  effort  over  the 
term  of  AFOSR  support  are  described  in  Section  2.  These  began 
as  a  fairly  narrow  interest  in  the  excitation  of  atmospheric 
gravity  waves  by  unstable  shear  layers,  but  rapidly  expanded 
to  encompass  many  aspects  of  gravity  wave  propagation, 
dissipation,  and  atmospheric  effects.  Of  particular  value  was 
the  application  of  a  theoretical  understanding  of  gravity  wave 


processes  to  the  interpretation  of  observational  data.  These 
findings,  together  with  the  results  of  analytic  and  numerical 
studies  performed  under  this  grant,  are  reviewed  in  Section  3. 
Research  under  this  grant  led  to  a  substantial  number  of 
publications  and  presentations.  A  cumulative  listing  of  these 
is  provided  in  Section  4.  Finally,  Section  5  lists  other 
peisonnel  participating  in  this  research  effort.  References 
are  provided  in  Section  6. 


2.  RESEARCH  OBJECTIVES 

The  initial  objective  of  research  under  t..is  AFOSR  grant 
was  to  understand  the  mechanisms  that  could  permit  the  rapid 
excitation  of  atmospheric  gravity  waves  by  unstable  shear 
layers.  This  excitation  is  thought  to  be  one  of  the  major 
sources  of  gravity  wave  motions  in  the  atmosphere  and,  as 
such,  must  be  better  understood  if  we  are  to  anticipate  the 
consequences  of  such  motions  on  the  large-scale  circulation 
and  structure  of  the  atmosphere.  Because  previous  linear 
studies  had  found  the  excitation  of  propagating  wave  motions 
to  be  constrained  by  very  small  growth  rates,  we  elected  to 
examine  the  effects  of  nonlinearity  on  this  process, 
specifically  the  excitation  of  large-amplitude  wave  motions 
via  the  nonlinefr  interaction  of  rapidly  growing  Kelvin- 
Helmholtz  instabilities. 
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A  second  objective  of  our  study  was  to  understand, 
through  numerical  modeling,  seme  of  the  consequences  of 
gravity  wave  saturation  in  the  middle  atmosphere.  These 
included  the  effects  of  local  convective  adjustment  of  the 
wave  field,  the  amplitude  limits  imposed  by  convective 
adjustment  for  single  waves  and  superpositions  of  waves,  the 
self  -  acceleration  of  gravity  wave  motions  due  to  induced  mean 
flow  changes,  the  excitation  of  other  motions  as  a  consequence 
of  local  adjustment,  and  the  relative  importance  of  nonlinear 
wave-wave  interactions  and  convective  adjustment  in  limiting 
wave  amplitudes.  The  goal  here  was  to  anticipate  what  wave 
processes  ought  to  be  most  significant  in  determining  the 
gravity  wave  spectra  and  characteristics  seen  in  atmospheric 
observations.  A  related  objective  addressed  the  effects  of  a 
localized  turbulent  diffusion  on  the  dissipation  of  wave 
motions  and  on  the  transport  of  heat  and  constituents. 

The  remaining  objectives  of  this  research  effort  dealt 
with  atmospheric  observations  and  their  interpretation.  A 
primary  focus  was  to  understand  the  very  different  echoes 
observed  by  the  Poker  Flat  MST  radar  in  the  mesosphere  during 
summer  and  winter.  The  radar  velocity  and  associated 
temperature  data  suggested  that  there  may  be  very  different 
mechanisms  at  work  in  the  two  cases.  A  related  goal  was  to 
identify  the  instabilities  responsible  for  wave  fi^.ld 
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saturation  and  the  turbulence  and  wave  amplitude  limits  they 
produce.  Also  of  interest  were  the  character  of  the  wave 
spectrum,  an  identification  of  the  motions  primarily 
responsible  for  energy  and  momentum  transports,  the 
consequences  of  local  turbulence  production  and  diffusion,  and 
the  relative  contributions  of  two-dimensional  (2-D)  turbulence 
and  gravity  waves  to  the  motion  spectrum  in  the  middle 
atmosphere . 


3.  SUMMARY  OF  RESEARCH  ACHIEVEMENTS 

Our  initial  research  objective  was  met  by  performing  a 
series  of  numerical  simulations  to  examine  the  excitation  of 
radiating  gravity  waves  both  via  the  linear  instability  of  a 
dynamically  unstable  shear  flow  and  through  the  interaction  of 
two  rapidly  growing  Kelvin - Helmhol tz  ( KH )  instabilities.  The 
former  simulations  revealed  that  radiating  waves  grow  in  a 
manner  consistent  with  linear  stability  analysis  until  the 
motions  reach  amplitudes  sufficient  to  cause  significant 
modifications  of  the  mean  flow.  Thereafter  the  wave  motions 
evolve  in  a  manner  that  is  dependent  on  the  wave  structure, 
with  motions  that  are  propagating  above  and  below  the  source 
shear  continuing  to  produce  significant  mean  flow  changes  at 
large  dis*-a.jces  from  the  shear  and  motions  that  are  evanescent 
above  or  *-alow  decaying  rapidly  due  to  a  continuing  erosion  of 
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the  shear  strength.  In  the  presence  of  nonlinearly  interacting 
KH  instabilities,  wave  motions  very  similar  to  the  linearly 
unstable  radiating  waves  were  found  to  be  excited  very 
efficiently,  to  achieve  large  amplitudes,  and  to  be  very 
transient  events  due  to  the  nature  of  the  excitation.  Because 
the  shear  geometry  in  this  study  allowed  for  the  presence  of 
radiating  unstable  modes,  however,  it  was  not  possible  to 
determine  whether  the  excitation  was  a  resonant  interaction  or 
simply  a  consequence  of  strong  nonlinearity  among  the  KH 
modes.  This  work  was  described  in  detail  by  Fritts  (1982a). 

A  subsequent  study  (Fritts,  1984a)  addressed  this  issue 
and  determined  that  the  latter  was  the  case,  with  the 
structure  of  the  radiating  modes  determined  not  by  linear 
instability  conditions,  but  by  the  nonlinearity  itself.  Also 
addressed  here  was  the  most  efficient  mechanism  for  such 
excitation.  That  found  to  be  most  important  was  the 
interaction  of  two  KH  modes  at  similar  horizontal  wavenumbers 
to  excite  a  radiating  wave  at  a  much  smaller  wavenumber 
(larger  scale).  A  second  mechanism  (vortex  pairing)  proposed 
by  other  authors  and  observed  in  unstratified  laboratory 
facilities  was  found  not  to  be  effective  in  the  stratified 
atmosphere . 

The  second  thrust  of  our  numerical  studies  vac  toward 
understanding  gravity  wave-mean  flow  interactions  ^.'ising  in 


response  to  wave  transience  and  dissipation  (including 
saturation).  A  first  study  addressed  one  evolution  of  a 
gravity  wave  packet  in  a  shear  flow  and  the  conditions  leading 
to  wave  instability  and  breakdown  (Fritts,  1982b).  It  was 
found  that  such  motions  can  be  strongly  influenced  by  a  time- 
dependent  mean  flow,  either  accelerating  or  retarding  the 
occurrence  of  instability,  and  that  when  such  instability 
occurs  for  high-frequency  motions  the  likely  consequence  is 
the  convective  instability  of  the  wave  field. 

Subsequently,  a  series  of  papers  (Dunkerton  and  Fritts, 
1984;  Fritts  and  Dunkerton,  1984;  Fritts,  1985)  examined  more 
specific  aspects  of  the  saturation  of  gravity  waves  in  the 
middle  atmosphere.  These  studies  utilized  a  convective 
adjustment  scheme  to  achieve  a  relaxation  of  the  unstable 
portions  of  an  evolving  wave  field  to  simulate  the  effects  of 
local  turbulence  generation.  The  significant  findings  here 
included  1)  an  amplitude  limit  largely  consistent  with  that 
expected  from  linear  saturation  theory  for  monochromatic  wave 
motions,  2)  a  substantial  acceleration  of  the  wave  phase 
speed,  denoted  se 1 f  -  accele r a t i on ,  due  to  its  residence  in  a 
region  experiencing  significant  mean  flow  accelerations  due  to 
wave  transience,  3)  a  saturated  amplitude  less  than  the 
monochromatic  limit  for  a  superposition  of  wave  motions,  4) 


the  excitation  of  additional  wave  motions  at  smaller  scales 


arising  as  a  consequence  of  the  local  adjustment  of  the  wave 
field,  5)  and  the  observation  that  nonlinearity  among  the 


saturating  wave  motions  is  not  sufficient  to  prevent  the 
occurrence  of  convectively  unstable  regions  within  the  wave 
field.  The  last  point  is  particularly  important  as  it  suggests 
that  gravity  wave  saturation  is  largely  a  linear  process  and 
does  not  require  nonlinearity  within  the  wave  field.  Point  2) 
implies  that  the  spectrum  of  gravity  wave  phase  speeds  in  the 
mesosphere  may  be  very  different  from  that  in  the  lower 
atmosphere,  and  points  1)  and  3)  are  consistent  with  recent 
observations  of  vertical  wavenumber  spectra  of  atmospheric 
motions.  Thus,  our  modeling  studies  have  given  us  some 
important  insights  into  the  mechanisms  responsible  for  the 
observed  motion  spectra  in  the  middle  atmosphere. 

In  a  closely  related  analytic  study,  Fritts  and  Dunkerton 
(1985)  examined  the  consequences  of  a  localized  turbulent 
diffusion  on  the  transports  of  heat  and  constituents  and 
concluded  that  such  an  effect  could  provide  for  significant 
wave  dissipation  while  dramatically  reducing  the  effective 
eddy  diffusion  acting  on  the  mean  thermal  and  constituent 
gradients.  It  is  believed  that  this  may  account  for  the  very 
different  observed  (Hocking,  1985)  and  required  (Strobel  et 
al.,  1986)  values  of  eddy  diffusion  in  the  mesosphere.  A 


second  analytic  study  completed  recently  ( Oy  et  al .  ,  1986) 


has  addressed  the  Stckes  drift  due  tc  vertically  prcpacating 
gravity  waves  and  has  shown  tnat  such  an  effect  can  account 
for  the  large  and  reversed  vertical  mean  motions  observed  by 

the  Poker  Flat  radar  relative  to  that  needed  to  balance  the 
meridional  circulation. 

More  recent  studies  under  this  AFOSR  grant  have  focussed 
on  the  interpretation  of  middle  atmosphere  motions  observed 
with  the  Poker  Flat  radar  (ie.,  Balsley  et  al . ,  1983).  Our 
initial  focus  was  on  the  large  asymmetry  between  the  summer 
and  winter  echoes.  We  now  believe  that  this  difference  is  due 
to  the  very  large  change  in  the  thermal  structure  of  the 
mesopause  region  with  season,  causing  winter  echoes  to  be 
characteristic  of  gravity  waves  saturating  largely  via 
convective  instability,  consistent  with  the  temperature  data 
of  Theon  et  al.  (1967).  Summer  echoes,  on  the  other  hand,  are 
concentrated  near  the  mesopause.  La r ge - ampl i tude  wave  motions 
are  less  likely  to  be  unstable  below  this  level,  but  are 
driven  to  saturation  amplitudes  by  the  rapid  increase  in 
stratification  at  this  height,  accounting  for  the  localized 
echoes  and  the  strong  increase  in  reflectivity  relative  to  the 
winter  condition. 

Also  examined  in  detail  was  data  collected  during  the 
STATE  experiment  (Fritts  et  al.,  1986).  This  study  used  both 


radar  and  rocket  velocity  and  temperature  data  to  define  the 


temporal  and  spatial  variability  of  the  m.tion  field  and 
address  the  mechanisms  responsible  for  turbulence  production 
within  the  wave  field.  It  was  found  that  the  motion  field  was 
dominated  by  a  large-scale  wave  motion  of  long  period  < ^  7  hr) 
that  strongly  controlled  the  intensity  and  location  of 
turbulence.  This  provided  strong  evidence  that  the  motion 
field  was  unstable,  with  a  wave  amplitude  near  that  required 
for  dynamical  instability  (more  likely  than  convective 
instability  for  low-frequency  motions)  and  with  turbulence 
occurring  at  that  location  in  the  wave  field  where  linear 
theory  predicted  instability.  Also  significant  here  was  the 
observation  that  small-scale  motions  similar  to  those  found  in 
the  numerical  study  by  Fritts  (1985)  were  present  in  those 
portions  of  the  large-scale  wave  field  expected  to  be  most 
unstable  and  thus  most  likely  sources  for  other  motions. 

Finally,  three  studies  addressing  the  wavenumber  and 
frequency  spectra  of  gravity  wave  motions  have  recently  been 
completed.  The  first  by  Smith  et  al.  (1985)  found  the  vertical 
and  oblique  wavenumber  spectra  to  be  consistent  with  a  gravity 
wave  model  due  to  VanZandt  (1985),  suggesting  that  gravity 
waves  are  a  much  more  significant  component  of  the  atmospheric 
motion  spectrum  than  is  2-D  turbulence.  A  second  study 
(Smith  et  al.,  1986)  showed  that  observed  vertical  wavenumber 
spectra  can  be  described  by  a  universal  saturation  spectrum 


with  a  saturation  power  consistent  with  linear  instability 
theory  and  the  wave  superposition  effects  noted  by  Fritts 
(1985).  This  is  potentially  very  significant  as  it  represents 
a  considerable  generalization  of  the  previous  monochromati c 
linear  saturation  theory.  The  final  spectral  study  by  Fritts 
and  VanZandt  (1986)  addressed  the  effects  of  Doppler  shifting 
of  gravity  wave  motions  and  showed  that  such  effects  can 
account  for  some  of  the  departures  of  observed  frequency 
spectra  from  those  predicted  in  the  absence  of  Doppler 
shifting  effects. 

Because  of  our  involvement  in  both  theoretical  and 
observational  studies  of  middle  atmosphere  gravity  waves,  we 
have  also  been  in  a  position  to  assess  and  review  developments 
in  the  field.  This  has  led  to  three  review  papers  dealing  with 
research  status  and  direction  in  this  field  which  hopefully 
have  had  a  positive  impact  (Fritts,  1984b;  Fritts  et  al., 
1984;  Fritts  and  Rastogi ,  1985). 


4.  PUBLICATIONS  AND  PRESENTATIONS 

During  this  four-year  research  effort,  AFOSR  has 
contributed  full  or  partial  support  for  25  presentations  at 
national  and  international  scientific  meetings,  including 
seven  invited  papers,  and  28  publications  that  have  appeared, 
are  in  press,  or  have  been  submitted  to  scientific  journals  or 
conference  proceedings.  The  presentations  include: 
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"Evidence  for  seasonally-dependent  generation  mechanisms,"  B. 

B.  Balsley,  W.  L.  Ecklund,  and  D.  C.  Fritts,  URSI  Int'l. 
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